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The determination of Os from hadronic r decays is impeded by the fact that two choices for the renormalisa- 
tion group resummation, namely fixed-order (FOPT) and contour-improved perturbation theory (CIPT), yield 
systematically differing results. On the basis of a model for higher-order terms in the perturbative series, which 
incorporates well-known structure from renormalons, it is found that FOPT smoothly approaches the Borel sum 
for the T hadronic width, while CIPT is unable to account for the resummed series. An example for the behaviour 
of QCD spectral function moments, displaying a similar behaviour, is presented as well. 
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1. INTRODUCTION 

Hadronic decays of the r lepton provide an ex- 
cellent playground for the study of QCD at low 
energies. Its mass of Mr « 1.8 GeV is in an en- 
ergy region where perturbation theory is still ap- 
plicable, but also non-perturbative effects come 
into play and need to be included. These may 
arise from vacuum-condensate terms in the frame- 
work of the operator product expansion (OPE) 
and from so-called duality-violating contributions 
close to the physical, Minkowskian energy axis. 

In the seminal article [T] , the strategy for a pre- 
cise determination of the QCD coupling from 
the total T hadronic width 



T[t — > hadrons /^t (7 ) ] 



T[r- 



■(7)] 



3.640(10), (1) 



was developed, while in the subsequent years mo- 
ments of spectral r-decay distributions were in- 
corporated into the analyses as well |2I3I4| . 

The analytical computation of the perturbative 
order correction [S] has recently revived the 
interest in analyses from hadronic t decays 
which after evolution to the Z boson mass scale 
resulted in the following determinations: 



a.(Af|) = 



0.1202 (6)oxp(18)th 

0.1212 (5)exp(9)th m 

0.1180 (4)exp(7)th m 

0.1187 (6)exp(15)th i 



(2) 



The dispersion in these results dominantly origi- 
nates from different treatments of the renormali- 
sation group (RG) resummation of the perturba- 
tive series, namely fixed-order perturbation the- 
ory (FOPT), or contour-improved perturbation 
theory (CIPT) |9|10j . being systematically larger 
than the last included term in the expansion, 
which often in asymptotic series provides an esti- 
mate of the uncertainty due to higher-order terms 
not included in the partial sum. This points to 
the necessity of investigating the influence of dif- 
ferent RG resummations in more detail, which 
can for example be performed within exactly re- 
summable models for the perturbative series. 

Most suitable for the as determination is the 
T decay rate into light u and d quarks Rt,v/a 
via a vector or axialvector current, since in this 
case power corrections are especially suppressed. 
Theoretically, Rt^v/a takes the form [I] 



IiT,V/A 



5, 



EwlKdP [l + <5(o) 



EW 



E e2 

D>2 



V/A 



(3) 



1.0198(6) [n] and 5'^ 



EW 



where Sew 

0.0010(10) [12] are electroweak corrections, 5*^°^ 
comprises the perturbative QCD correction, and 
the ^^^y/^ denote quark mass and higher D- 
dimensional operator corrections which arise in 
the framework of the OPE. 
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2. PERTURBATIVE CORRECTION S^°^ 

Below, only the purely perturbative correction 
shall be considered, which gives the dominant 
contribution to Rry/A- In FOPT it takes the 
general form 

oo n 
ra=l fc=l 

where a(/i^) = = as(/i)/7r, and c„^fc are the 
coefficients which appear in the perturbative ex- 
pansion of the vector correlation function, 

nv(^)--Y^E<E^-^i-' 7^ ■ (5) 

n=0 fc=0 \ / 

At each perturbative order, the coefficients c„_i 
can be considered independent, while all other 
Cn,k with k>2 are calculable from the RG equa- 
tion. Further details can for example be found 
in ref. [?]■ Finally, the J; are contour integrals in 
the complex s-plane, which are defined by 

1 /■ dr 

Ji^ TT- f — (l-x)3(l + x)ln'(-x). (6) 

\x\ = l 

The first three, being required up to 0{al), take 
the numerical values 

7 _ 1 J _ 19 T _ 265 1 _2 /y^i 

At order a" FOPT contains unsummed loga- 
rithms of order ln'(— x) ^ tt' with I < n related 
to the contour integrals Ji . CIPT sums these log- 
arithms, which yields 

oo 
n=l 

in terms of the contour integrals J^^{M^) over the 
running coupling, defined as: 

J^iM^) ^ ^ / ^ (1-^)' (l+x)a"(-M» . 

(9) 

In contrast to FOPT, for CIPT each order n 
just depends on the corresponding coefficient Cn,i- 



Thus, all contributions proportional to the coeffi- 
cient c„^i which in FOPT appear at all perturba- 
tive orders equal or greater than n are resummed 
into a single term. 

Numerically, the two approaches lead to sig- 
nificant differences. Employing the recent aver- 
age UsiMz) = 0.1184 nil, leading to as(M^) = 
0.3186, in eqs. (g]) and ©, one finds 

4o = 0.1959 (0.2022) , (10) 

S^ci = 0.1814 (0.1847), (11) 

where the first number in both cases employs the 
known coefficients up to ©(a^) [5] and the num- 
bers in brackets include an estimate of the C(q;^) 
term with cs^i « 283 7^. Inspecting the individ- 
ual contributions from each order, up to 0{a^) 
the CIPT series appears to be better convergent. 
However, around the seventh order, the contour 
integrals J^(M^) change sign and thus at this or- 
der the contributions are bound to become small. 
Therefore, the faster approach to the minimal 
term does not necessarily imply that CIPT gives 
the closer approach to the true result for the re- 
summed series. 

3. A PHYSICAL MODEL 

To investigate whether FOPT or CIPT results 
in a better approximation to S^^\ one requires a 
physically motivated model for its series. Such a 
model was constructed in ref. [7] and is based on 
the Borel transform of the Adler function Dv{s): 

Dv{s) EE -S^ny(s) EE ^ [1 + D{s)] . (12) 

In the following discussion it is slightly more con- 
venient to utilise the related function D{s). Its 
Borel transform B[D]{t) is defined by the relation 

oo 

D{a) = J dte-*/°'B[D]{t). (13) 



The integral D{a), if it exists, gives the Borel sum 
of the original divergent series. It was found that 
the Borcl-transformed Adler function B[D]{t) ob- 
tains infrared (IR) and ultraviolet (UV) renor- 
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malon poles at positive and negative integer val- 
ues of the variable u = 9f/(47r), respectively 
[14115116) . (With the exception of it = 1.) 

Apart from very low orders, where a dominance 
of renormalon poles close to m = has not yet set 
in, intermediate orders should be dominated by 
the leading IR renormalon poles, while the lead- 
ing UV renormalon, being closest to u = 0, dic- 
tates the large-order behaviour of the perturba- 
tive expansion. Assuming that only the first two 
orders are not yet dominated by the lowest IR 
renormalons, one is led to the ansatz 

B[D]{u) = B[Dr]{u) + Bi&J^'Ku) + 

+ + (14) 

which includes one UV renormalon at u = —1, the 
two leading IR renormalons at u = 2 and u — 3, 
as well as polynomial terms for the two lowest per- 
turbative orders. Explicit expressions for the UV 
and IR renormalon pole terms B[Dp^]{u) and 
B[D^]{u) can be found in section 5 of ref. 

Apart from the residues d^^ and cJ^^, the full 
structure of the renormalon pole terms is dic- 
tated by the OPE and the RG. Therefore, the 
model p4|) depends on five parameters, the three 
residua d^^, ^2^^ and as well as the two poly- 
nomial parameters dg and d^'^ . These parame- 
ters can be fixed by matching to the perturbative 
expansion of D{s) up to 0{al). Thereby, also the 
estimate for C5.1 is used. The parameters of the 
model (fT4|) are then found to be: 

d^^ = -1.56-10-2, 4^= 3.16, 4^= -13.5, 

4° = 0.781, 4° = 7.66- 10-^ (15) 

The fact that the parameter df'-' turns out to be 
small implies that the coefficient C2,i is already 
reasonably well described by the renormalon pole 
contribution, although it was not used to fix the 
residua. Therefore, one could set ° = and 
actually work with a model which only has four 
parameters. The predicted value cs^i — 280 then 
turns out very close to the estimate, which can be 
viewed as one test of the stability of the model. 

This is also corroborated in table 1, where 
the relative contributions of a certain renormalon 



C2a C34 C4^i €5,1 C6,i 

IR2 -77.8 82.4 100.4 135.9 97.5 

IR3 152.0 28.7 -10.0 -20.2 -13.3 

UVi 22.5 -11.2 9.7 -15.6 15.8 



Table 1 

Relative contributions (in %) of the different IR 
and UV renormalon poles to the Adler-function 
coefficients C2,i to cg,! for the Borel model (ITU) . 

pole to the coefficients C2.1 to cg,! is tabulated. 
The sum of the contributions to C2,i is close to 
100%, implying again that the polynomial term is 
already small. Then, from C3_i to cg,! the leading 
IR pole at M 2 is dominating the coefficients, 
before the leading UV pole at m = — 1 takes over 
at even higher orders. (For the central model (|14p 
this happens around the 10th order [7].) 
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Figure 1. Results for (5pQ (full circles) and i5^j 
(grey circles) at OLg^Mr) = 0.3186, employing the 
model dU]), as a function of the order n up to 
which the terms in the perturbative series have 
been summed. The straight line represents the 
result for the Borel sum of the series. 

The implications of the model for (5*^°^ in 
FOPT and CIPT is graphically represented in fig- 
ure [H The full circles denote the result for S^q 
and the grey circles the one for ^^j-*, as a func- 
tion of the order n up to which the perturbative 
series has been summed. The straight line cor- 
responds to the principal value Borel sum of the 
series, Jg^g = 0.2080, and the shaded band pro- 
vides an error estimate based on its imaginary 
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part divided by tt. The order at which the series 
have their smahest terms is indicated by the grey 
diamonds. As is obvious from figure [U FOPT 
displays the behaviour expected from an asymp- 
totic series: the terms decrease up to a certain 
order around which the closest approach to the 
resummed result is found, and for even higher or- 
ders, the divergent large-order behaviour of the 
series sets in. For GIFT, on the other hand, the 
asymptotic behaviour sets in earlier, and the se- 
ries is never able to come close to the Borel suml3 





C2,l 


C3,l 


C4,l 


C5,l 


C6,l 


IR3 


-743.3 


-140.5 


49.1 


98.9 


99.1 


IR4 


662.8 


244.2 


47.7 


6.3 


-7.2 


UVi 


7.5 


-3.7 


3.2 


-5.2 


8.1 



Table 2 

Relative contributions (in %) of the different IR 
and UV renormalon poles to the Adler-function 
coefficients C3 1 to ce 1 for the Borel model with 
df =0. ' 



As the behaviour of GIFT versus FOFT hinges 
on the contribution of the leading IR renormalon 
at M = 2, in principal also models can be con- 
structed for which GIFT provides a better ac- 
count of the Borel sum. These would generally 
be models where is much smaller than the 
value quoted in eq. ([15]). While such models can 
at present not be excluded, the pattern of the 
individual contributions appears more unnatural 
than in the main model (jl4p : the known Cn,i can 
only be reproduced when one allows for large can- 
cellations between the individual terms. Thus, 
the behaviour generally expected from the pres- 
ence of renormalon poles, namely dominance of 
leading IR poles at intermediate orders, would be 
lost. 

This is apparent from table [2l which is the ana- 
log of table [11 but for a model where is forced 
to be zero and an additional IR pole at u = 4 
is added, in order to be able to reproduce the 
known Adler function coefficients. On the one 

^The same conclusions had already been drawn in ref. 1171 
on the basis of the large- /Jq approximation for Ilv{s). 



hand for low coefficients there are huge cancella- 
tions between the IR renormalon poles and also 
C2_i is not well described at all. This entails that 
a large, additional polynomial term is required. 
Besides, it appears unnatural that the residue of 
the first IR renormalon pole at u = 2 is small, and 
still there should be a natural size contribution of 
the gluon condensate to the Adler function. 

A graphical account of the model with — 
is presented in figure [5| As anticipated, now 
GIFT provides a good description of the Borel 
sum, while FOFT is able to come reasonably close 
to it around its minimal term, but generally is 
rather badly behavedH Nonetheless, again, the 
behaviour observed in table [D and figure [5] ap- 
pears unnatural from the perspective of the struc- 
ture of the Borel transform of the Adler function 
and should be considered less likely than the be- 
haviour of the model ([T4| with the residues (|15p . 



0.22 - 




Perturbative order n 



Fi gure 2. Results for S^q (full circles) and 
(grey circles) at asiM^.) = 0.3186, employing the 
model ([Ti)) with o?2^ = and an additional IR 
pole at u = 4, as a function of the order n up to 
which the terms in the perturbative series have 
been summed. The straight line represents the 
result for the Borel sum of the series. 

In the standard as determinations from 
hadronic r decays [2|3| . besides the total decay 
rate also moments of the spectral decay distri- 
butions are employed. Historically, the so-called 

similar behaviour was found in ref. 18 in models where 
the higher-order Adler function coefficients were assumed 
to be small, in contrast to the expected asymptotic be- 
haviour of this series in QCD. 
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(fc, Z)-moments were used, for which a polynomial 
(1 — x)''x'' is multiplied to the kinematical weight 
function from phase-space. In the analyses [213] . 
moments with k = 1 and / = 0,1,2,3 were taken 
into account, such that also condensate contri- 
butions up to dimension-8 could be extracted in 
addition to as- Therefore, it is of interest to anal- 
yse the behaviour of these moments in models of 
higher orders of perturbation theory as well. 




Perturbative order n 



Figure 3. Results for (1,2) moments S^q (full 

circles) and S''^{^^^ (grey circles) at as{Mr) — 
0.3186, employing the model (fT4|) . as a function 
of the order n up to which the terms in the per- 
turbative series have been summed. The straight 
line represents the result for the Borel sum of the 
series. 

An example of such an analysis is shown in fig- 
ure |3] for the moment (1,2). It is observed that, 
like for 5'"^ CIPT is unable to provide a rea- 
sonable account of the full Borel sum. In con- 
trast, FOPT comes close to the resummed result 
for perturbative orders around the minimal term, 
though it is obvious that this particular moment 
displays a very bad behaviour of the asymptotic 
series, with only unsatisfactory convergence up 
to its minimal term. The example of the (1,2)- 
moment should motivate an exhaustive investi- 
gation of the moments employed in previous as 
analyses from hadronic r decays, which will be 
presented in the near future. 

4. CONCLUSIONS 

Models of higher orders of perturbation the- 
ory allow for the study of different resummation 



prescriptions in the computation of the total r 
hadronic width, as well as related moments of 
T decay spectral distributions. The most promi- 
nent methods are fixed-order perturbation theory 
(FOPT) and the so-called contour-improved per- 
turbation theory (CIPT), which performs a par- 
tial resummation of running effects of the QCD 
coupling as in the integration along the complex 
contour in the s-plane. 

A physically motivated model for the higher- 
order behaviour of the Adler function was pre- 
sented in ref. [7 , and is given in eq. The 
model was based on the general structure of the 
Borel transform of the Adler function and the 
renormalisation group equation. Furthermore, 
knowledge on the first four analytically available 
coefficients ci^i to c^^i, as well as an estimate for 
the fifth coefficient cs^i, were incorporated. 

Results for ^'^^ in the main model ([T4l) were dis- 
played in figure [U and it is observed that while 
FOPT provides a good account of the full Borel 
summation, CIPT is never able to come close to 
the resummed valued This general behaviour 
hinges on the size of the residue of the first IR 
renormalon pole at m = 2. In a model in which 
this residue is set to zero by hand, on the con- 
trary CIPT well describes the Borel sum, whereas 
FOPT, though approaching the resummed value 
around its minimal term, generally is rather badly 
behaved. Models in which ~ 0, however, are 
only able to reproduce the known Adler function 
coefficients through large cancellations between 
different IR renormalon contributions, which ap- 
pears unnatural. 

In the standard experimental extractions of as 
from hadronic r decays, also moments of the de- 
cay spectra are employed |2I3) . An example of 
the behaviour of such a moment was shown in 
figure [H The repeatedly unacceptable behaviour 
of CIPT in this case and the also unsatisfactory 
convergence of FOPT, suggest that such moments 
should be investigated systematically, before their 



^The general behaviour of the Borel model is also sup- 
ported by an independent approach where the perturba- 
tive series is conformally transformed into a series which 
displays better convergence properties than the original 
series in powers of as |19j . though in this case a modified 
CIPT better converges towards the full result. 
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usefulness in as determinations from hadronic r 
decays is corroborated. 

A final topic, not touched upon at all so far, 
are violations of quark-hadron duality [20121122] . 
Like the perturbative expansion in a^, also the 
OPE in inverse powers of s may be only asymp- 
totic, and thus exponentially suppressed terms 
could be relevant close to the Minkowskian axis 
where the bound states are situated. This was 
investigated in a model in [51] and the possible 
influence of duality violations in t decay spectra 
was studied in refs. [12] ■ Also here a more sys- 
tematic investigation seems to be in order which 
was initiated in [23j and will be continued in the 
future. 

Hadronic r decays have proven to be a very 
fruitful laboratory for the study of low-energy 
QCD and the extractions of fundamental QCD 
parameters like the coupling a^.. However, there 
remain unresolved theoretical issues which taint 
the precision of these determinations. Numeri- 
cally the two most relevant appear to be the re- 
summation of QCD running effects in the compu- 
tation of the T hadronic width and related decay 
moments as well as duality violations. As far as 
the former topic is concerned, a physically moti- 
vated model of the QCD Adler function favours 
the use of FOPT, since generally it is better able 
to represent the fully resummed series. 
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